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The primary target for the planned space-borne gravitational wave interferometers DECIGO /BBO 
is a primordial gravitational wave background (PGWB). However there exist astrophysical fore- 
grounds and among them, gravitational waves from neutron star (NS) binaries are the solid and 
strong component that must be identified and subtracted. In this paper, we discuss the geometry 
of detector configurations preferable for identifying the NS/NS binary signals. As a first step, we 
analytically estimate the minimum signal-to-noise ratios (SNRs) of the binaries for several static 
detector configurations that are characterized by adjustable geometrical parameters, and determine 
the optimal values for these parameters. Next we perform numerical simulations to take into ac- 
count the effect of detector motions, and find reasonable agreements with the analytical results. We 
show that, with the standard network formed by 4 units of triangle detectors, the proposed BBO 
sensitivity would be sufficient in receiving gravitational waves from all the NS/NS binaries at z < 5 
with SNRs higher than 25. We also discuss the minimum sensitivity of DECIGO required for the 
foreground identification. 



o 

u 

43 

Or 

6 



PACS numbers: Valid PACS appear here 



I. INTRODUCTION 



Currently, several ground-based gravitational wave interferometers (e.g. LIGO, VIRGO, GEO and TAMA) are 
operating or in installation/commissioning phases for upgrades. In the next decade, second generation detectors (e.g. 



advanced-LIGO, advanced- Virgo and LCGT) would be available. These have optimal sensitivities around 100-1000 

o : 

, Hz aiming to detect gravitational waves (GWs) from NS/NS inspirals, supernovae, etc. Furthermore, the possibilities 

ov 

of more powerful third generation detectors (e.g. ET) have been actively discussed with primordial gravitational wave 

t-H ■ 

background (PGWB) as one of their principle targets, and they would newly explore the lower frequency regime down 

to ~ 1 Hz. Meanwhile, NASA and ESA are planning to launch a space-borne gravitational wave interferometer called 

> ' IT 

the Laser Interferometer Space Antenna (LISA) [1] whose optimal band is 0.1 — 100 mHz. GWs from white dwarf 

(WD) binaries are guaranteed targets for LISA, but they will also form confusion noises that might mask other signals 

including PGWB (see Refs. 2|, y] for reviews on this source). 

The Deci-Hertz Interferometer Gravitational Wave Observatory (DECIGO) is a future plan of Japanese space 



mission for observing GWs around / ^ 0.1 — 10 Hz 



4J-|6|. In US, the Big Bang Observer (BBO) has been proposed 

for a proposal 



y, 



as a follow-on mission to LISA, and its optimal band is similar to that of DECIGO [7| (see also Ref. 
of an atomic interferometer). Since GWs from WD/ WD binaries have cut off frequencies at / < 0.2Hz [9|, the higher 
frequency part of DECIGO/BBO is free from their confusion noises. Therefore, the direct detection of the PGWB 
produced during inflation has been set as the primary goal of DECIGO and BBO. 



For the amplitude of PGWB generated by standard slow-roll inflation, the WMAP team has placed a conservative 
bound corresponding to f^cw J$ 10~ 14 around 1 Hz 



10j . where &Gw(f) represents the normalized energy density of 
GW background at a frequency /. However, with theoretical analysis based on the temperature anisotropics reported 
by WMAP, we can reasonably expect a level Slew ^ 10~ 15 for the inflation background at the DECIGO/BBO 
band. Since the amplitude of the target PGWB is expected to be small, the correlation analysis would be a powerful 
technique for its detection 
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12| . The sensitivity of this method is characterized by the so-called overlap reduction 
function that depends strongly on the geometry of the detector network. It becomes maximum for a pair of co-aligned 
detectors. With DECIGO and BBO, two units of aligned detectors in the star-of-David form are planned to be used 
for correlation analysis with their detection limit ficw ~ 10~ 15 -10~ 16 . In addition to these co-aligned ones, two other 
units might be launched, mainly to improve the angular resolution of astrophysical sources, by forming long-baseline 
(~1 AU) interferometers (see Fig. [I]). 

While the frequency regime of DECIGO and BBO would not be fundamentally limited by WD/WD binaries, NS/NS, 
NS/BH and BH/BH binaries (BH: black hole) are considered as their important astrophysical sources. Among them, 
the foreground GWs by NS/NS binaries are understood relatively well with their cosmological merger rate (10 5 /yr) 



estimated observationally 



13, 



14[. Therefore, for successful 



14|. Its total amplitude is expected to be I^gw ~ 10~ 12 
operations of DECIGO and BBO in opening deep windows of GWs around 1 Hz, it is essential to identify NS/NS 
(also NS/BH and BH/BH) binaries and remove their contributions from data streams of detectors by ^4-5 orders of 
magnitude in terms of ^gw- 



Cutler and Harms 14| have done a pioneering work on this topic. They set the threshold signal-to-noise ratios 
(SNRs) pthr for binary detection and self-consistently estimated the fraction of binaries whose SNRs are below jOthr- 
This corresponds to the fraction of residual NS /NS binaries that cannot be subtracted. According to their calculation, 
with the base design sensitivity of BBO, it is possible to subtract out NS/NS binary signals more than 5 orders of 
magnitude. However, if it is worse by a factor of 2, the prospect crucially depends on the value of pthr- If the 
sensitivity is worse by a factor of 4, then it seems difficult to clean the foreground noises down to the level needed 

n 

for the detection of the inflation background. Harms et al. [15| have performed numerical simulations of a projection 



14|. 



method for reducing the subtraction noises and confirmed the results obtained by Cutler and Harms 

Our basic aim in this paper is to specifically study how the performance of the binary identification depends on 
the geometry of the detector network of DECIGO or BBO. Here we should comment on some of the interesting 
aspects on the detector configuration. If there is no astrophysical foreground, in order to detect PGWB, it is most 
favored to put all the detectors co-aligned at the same location so that the correlation analysis would work for all the 



pairs of detectors. However, in order to identify NS/NS binaries and subtract them, detectors should have different 
orientations to cover up the whole sky. This is particularly relevant for a binary at a high redshift. Such a binary goes 
through the DECIGO/BBO band in a short amount of time (e.g. ~1 month for a NS/NS binary at z ~ 5), and we 
need multiple detectors with different orientations not to miss its detection. Thus, determining each configuration of 
DEC1GO and BBO is an important topic closely related to their primary scientific goal. Since the design parameters 
of DECIGO are currently under discussion jl 61 ] . another goal of this paper is to clarify its sensitivity required to 
sufficiently identify the NS/NS binaries. 

We consider several examples of possible detector networks by introducing geometrical parameters that should be 
adjusted to maximize the minimum value of the signal-to-noise ratio (SNR) of the binaries. The minimum value 
would be a helpful reference to discuss the prospect for the identification of all the foreground binaries, and we heavily 
rely on this reference. In addition to simple but workable analytical evaluations for the minimum SNRs with static 
detector networks, we include the effect of detector motions by performing Monte Carlo simulations with the code 



developed in Refs. 



17 



18(. Then we discuss the detector sensitivity required for identifying all the NS/NS binaries. 
Our analysis might also provide useful insights on the geometrical properties of ground-based detector network. 

In this paper, we only discuss the GW foreground made by NS/NS binaries, as a concrete example. We can expect 
that BH/NS or BH/BH binaries would cause less severe problems, because of their larger chirp masses (i.e. larger 
signals) and their merger rates presumed to be smaller than NS/NS binaries. 

This paper is organized as follows. In Sec. [TTJ we briefly describe GW observation with DECIGO and BBO. In 
Sec. IIII1 we discuss the basic properties of GW foreground of NS/NS binaries. In Sec. IIV1 we study the responses 
of detectors to incoming GW from a binary. In Sec. [V] we explain the dead angles for GW interferometers and 
analytically calculate the minimum SNRs for various static detector configurations. In Sec. IVI1 we first explain the 
set ups for our Monte Carlo simulations. Then, we provide the numerical results and their interpretations. In Sec. lVIll 
we summarize our work, suggesting the preferable configurations for DECIGO and BBO, and comment on possible 
future works. In Appendix we provide formulae that would be useful when relating angular parameters in a 
spatially fixed frame with the ones in a rotating frame attached to a detector. 
We take the unit G = c = 1 throughout this paper. 
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FIG. 1: The default configuration of BBO or DECIGO with az = 120°. Totally four units of triangle-like detectors will be 
operated. Two of them are nearly co-aligned to form a star-of-David constellation for the correlation analysis. Two outrigger 
ones are used to improve localization of individual astrophysical sources. 

II. DECIGO/BBO 
A. Two effective interferometers 

In its original proposal, BBO consists of four units of detectors. Each unit is composed of three drag-free spacecrafts 
to form a nearly regular triangle. Its default configuration is shown in Fig. [1] with 03 = 120° 7]. These four units 
are tilted 60° inwards relative to the ecliptic plane to keep its arm lengths nearly constant (as for LISA), and move 
around the sun with the orbital period of 1 yr. A similar configuration is supposed to be adopted for DECIGO. The 
main difference between these two is the basic design of detectors. While DECIGO is planned to be Fabry-Perot type 
interferometers, BBO will be transponder-type interferometers like LISA. 

Here, we briefly sketch a rough idea to generate independent data streams from each triangle-like unit of BBO (see 
e.g. Ref. [19| for a detailed analysis with LISA). First we make three interferometers a, b and c defined at the three 
vertices by using their adjacent two arms (with opening angles 60°, see Fig. [5]) symmetrically. Since they share same 
arms of interferometers, noises in the data from these three interferometers have correlations. Here, we assume that 
the covariance matrix of these noises has a simple symmetric structure. Namely, its diagonal elements are the same, 
and the off-diagonal elements have an identical value. Then, in order to obtain the orthogonal data streams, we take 
the linear combinations I and II as follows; 

I= a + b_2c n= a_b 
y/6 y/2 



These essentially correspond to the A and E modes of the TDI data streams in Ref. 



The T mode is not important 
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FIG. 2: One unit of BBO (DECIGO) detector and the orientations of two effectively L-shaped interferometers I and II defined 
in Eq. (1}. 

for our study. Due to the symmetry of the noise matrix of the original data a, b and c, the noises of the combined data 
I and II have an identical spectrum with no covariance between them 



1JH22| 



In the long- wave limit, responses of 



I and II to GWs can be effectively regarded as those of two L-shaped interferometers whose orientations are shown in 



Fig. m 



21 



23|. The interferometer II is obtained by 45° rotation of the interferometer I. Note that the I-II pair forms 
an orthogonal basis for L-shaped interferometers on the detector plane. In this paper, we apply the above arguments 
also for each (triangle-like) unit of DECIGO, assuming that the noise spectra of its two effective interferometers are 
identical and uncorrelated. With the four units of triangles as shown in Fig. [1] the total number of the effective 
(L-shaped) interferometers becomes eight. 

B. DECIGO/BBO noise spectrum 

In this subsection, we provide the noise spectrum of DECIGO and BBO for each effective interferometer mentioned 
in the previous subsection. The data stream s a (t) of an interferometer a (e.g. a —I or II) can be decomposed into 
GW signal h a (t) and the detector noise n a (t) as 



s a (t) = h a (t) + n a (t). 



(2) 



In Sec. IV, the concrete expressions of the former h a (t) are presented for individual chirping binaries. Assuming that 
the noise n a (t) is statistically stationary, we take its Fourier transform as 



a (.f) = / dte^n^t), 



(3) 



and define its spectrum Sh(f) by 



(n a (f)nMT) = ^(f ~ f)S h (f). 



(4) 



Here (•) represents an ensemble average. Since we only use the effectively L-shaped interferometers whose noise 
spectra are identical for given missions (DECIGO and BBO), we omitted the label a for the spectrum. In this paper 
we conservatively deal with an appropriate frequency regime / > 0.2Hz, considering the potential effects of WD/WD 
binaries 

For its reference design parameters 



-.DECIGO 



(/) = 6.53 x 10 



-49 



1 



371 ]. the noise spectrum of DECIGO is fitted as 

/ v 4 i 



+ 4.45 x 10 



-51 



1 Hz 



1 



+ 4.94 x 10 



-52 



1 Hz 



Hz" 1 , (5) 



with f p = 7.36Hz. The first, second and third terms represent the shot noise, the radiation pressure noise and the 
acceleration noise, respectively. We have checked that, in the frequency regime 10~ 2 -10Hz, this formula reproduces 



the sensitivity curve in Fig. 2 of Ref. 5] with an error less than 1.41% [36 

n 

For BBO with its reference parameters [39j, we fit its noise curve shown in [14[ and use the following (non sky- 
averaged) expression; 



^ B °(/) = 2.00x 10 



-49 



1 Hz 



4.58 x 10 



-49 



1.26 x 10- Hz- 



(6) 



that is almost equivalent to the spectrum used in Cutler and Holz [24| (up to a factor 5 corresponding to sky- 
averaging). In Fig[3l we present these noise spectra for each L-shaped interferometer. Here, in order to compare with 
the characteristic amplitude of a NS /NS binary, we show the sky averaged forms of the spectra which are a factor of 
5 larger than the original ones expressed in Eqs. (0 and (|6]). As shown in Fig. [3j the noise spectrum of DECIGO is 
~ 3 times larger than that of BBO. We mainly use the spectral shape of DECIGO for our numerical evaluation below 
and later we discuss its required level for sufficient NS/NS cleaning. We also show that, in appropriately normalized 
forms, our results for DECIGO are almost the same as those for BBO. Therefore, our analyses can be easily applied 
to BBO as well. 
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FIG. 3: The solid curves show the noise spectra for DECIGO (thick red) and BBO (thin blue). Note that these noise curves are 
the sky-averaged ones (a factor of \/5 larger than the original ones) . The (black) thin dotted line represents the expected total 
amplitude of NS/NS foreground and the (purple) thick dotted line represents the primordial GW background corresponding to 
JIgw = 10 -16 . The (black) dashed line at / = 0.2Hz indicates the upper frequency cutoff of WD/WD confusion noises. 



III. GRAVITATIONAL WAVE BACKGROUND FROM NS/NS BINARIES 



Magnitude of a stochastic GW is characterized by its fractional energy density JIgw per logarithmic frequency 
interval as 



SI 



GW 



1 dp GW 
Pc dlnf 



(7) 



where pew is the energy density of the GWs and p c = is the critical energy density of the universe with the current 
Hubble parameter Hq = 70/i7o km/s/Mpc with hn = 1 for our fiducial value 



cosmological NS/NS binaries, we apply the convenient formula given by Phinney 26j as 

n(z) 



*25j. For tin- total GW foreground by 



Sir 5 / 3 1 



9 Hi 



dz- 



(8) 



/ {l + z)*/*H(z)' 

Here M. is the chirp mass defined as M. = (mini2) 3, ' 5 (™i + m^ 1 ^ with the two masses m\ and ni2 of binaries. In 
this paper, we assume mi = m 2 = 1.4M Q (M Q = 2.0 x 10 33 g) for NS/NS binaries. The function n{z) is the NS/NS 
merger rate per proper time per comoving volume at redshift z, and the Hubble parameter H(z) is given by 



H(z) = H Vn m (l + z) 3 + n A 



(9) 



for the ACDM cosmology assumed in this paper with the cosmological parameters Jl m = 0.3 and J1a = 0.7. We 
re-express the rate h{z) in terms of the current merger rate no and the redshift dependence s(z) as n(z) = uq x s(z). 
In this paper, we use tiq = 10~ 7 Mpc~ 3 yr _1 as a fiducial value. For the redshift dependence s(z), we adopt the 



following piecewise linear fit 



14 



27] based on Ref. 



28|; 



1 + 2z (z < 1) 

*(*)=< |(5 -«) (1<«<5) (10) 
(z>5). 



Then, the magnitude f2 G w is numerically evaluated as 



^gw(/) - 1-73 x io- 12 ^ 3 ( _^_y /a ( j_y 3 c — «o — ) . ,i] , 

70 \1.22MqJ \mzj \lO-mpc- 3 yr- l J 



The total (sky-averaged) GW foreground spectrum Sh and the normalized energy density ficw have the relation 



14] 



and, with Eqs. (JTTJ and (|12p. we obtain 



- -.TV^gw, (12) 

7T 



5/6 / . \ 1/2 / , v -7/6 



,/^(/) = 1.17xlO- 2 Vo 1/2 (— —J ( ^ J ( — J Hz" 1 / 2 . (13) 



For a PGWB, we have 



n GW (f)\ 1/2 f f y 3/ \ _! 



5 PGWB (/) = 8 . 85xl0 -27^ o ^^Uij Hz -l/2_ (14) 

These relations are also shown in Fig. El As mentioned in the introduction, to detect a PGWB as small as Oqw ~ 10 -16 
around / ~ 1Hz, we need to clean the NS/NS foreground and reduce its residual by 5 orders of magnitude in terms 
of tt G w- 

IV. GW FROM BINARY AND DETECTOR RESPONSE 

In order to subtract the GW foreground formed by cosmological NS/NS binaries, it is essential to individually 
identify them in the data streams of detectors. Here, we summarize the basic formulae required for evaluating the 
SNRs of these binaries. 

First, we provide expressions for the GWs coming from a NS/NS binary. Since we start our calculation from 0.2 
Hz where eccentricity of the binary is expected to be small, we can reasonably assume a circular orbit for studying 
detectability of the binary. From the quadrupole formula of GWs, the waveforms of + and x modes in the principal 
polarization coordinate are given by 

h + (t) = A + cos (15) 
h x (t) = A x sm4>(t), (16) 
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FIG. 4: We use two types of coordinates: (i) a barycentric frame (x, y, z) tied to the ecliptic and centered in the solar system 
barycenter, (ii) an detector frame (x, y, z) attached to the detector as in Fig. [21 



with the phase function <p(t) and the amplitudes 

A+ = 

A-* = 



rD L 
4mi?ri2 



rD L 



(! + (£• AT) 2 ), 
L ■ N). 



(17) 
(18) 



Here r is the orbital separation of the binary, Dl is the luminosity distance to the source, L is the unit vector parallel 
to the orbital angular momentum and N is the unit vector pointing towards the source from the observer (see Fig. 2]). 
Next, we deal with responses of GW detectors. Following Ref. 21| . we introduce two Cartesian reference frames: 



(i) a barycentric frame (x, y, z) tied to the ecliptic and centered in the solar system barycenter, with z (unit vector 
in z direction) normal to the ecliptic, (ii) an detector frame (x, y, z) attached to the detector, with the direction z 
normal to the detector plane (see Figs. [2] and S]). 

While each unit of DECIGO or BBO changes its position and orientation as a function of time, we first discuss its 
response by fixing the geometry of the source-detector system. As already discussed in Sec. Ill Bt the triangle unit can 
be effectively regarded as the two L-shaped interferometers I and II. Each output of these interferometers is written 
as 



h a {t) = A + F+(9 S , <p s ,ip 3 ) cos cj>{t) +A X F* (0 S , 0s, V>s) sin 0(t) 



(19) 
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with the label for interferometers a=I, II. The beam pattern functions for the interferometer I are given by 



F+(9 s ,cj)s, if>s) = -(1 + cos 2 S ) cos(20 s ) cos(2Vs) - cos(8 s ) sin(20 s ) sin(2V>s), 
Fj x (0 S , fa, Vfe) = 2^ + cos2 ^s) cos ( 2 ^s) sin(2Vs) + cos(0 s ) sin(20 s ) cos(2?/>s), 



and those for the interferometer II are expressed as 



F*(0 s ,<£ s ,Vs) = F*(fiB,<h-*/*,ik)- 



(20) 
(21) 

(22) 
(23) 



Here, the two angles (0s,4>s) represent the direction of the source in the detector frame and is the polarization 
angle given as 

L-z-(L- N)(z ■ N) 



tan^s 



N-(Lxz) 



(24) 



This angle characterizes the orientation of the vector L projected on the transverse plane (see Fig. 1 of Ref. 291 ]) 
By using the stationary phase approximation, we obtain the Fourier transform of the signal as 



h a (f) = Af 



= A f- 7 /6 P **(/) 



5 

^ ^pol,a 



-i(,tPpal,a+<P d) 



where the amplitude A is given by the redshifted chirp mass M z = (1 + z)M as 



A = 



1 M 



5/6 



307r 2 /3 D L 



The polarization amplitude A po i ta is defined by 



A po l,a =J(l + (L- iV)2)2(F+)2 + 4(i • AT)2(F Q X )2. 



(25) 



(26) 



(27) 



In Eq. (|25[) , is the phase in frequency domain, ip po i,a is the polarization phase, and ifo is the Doppler phase 

which denotes the difference between the phase of the wavefront at the detector and the one at the solar system 
barycenter. Since we are only interested in the signal to noise ratio (SNR) of binary GWs, we do not need explicit 
forms for the phases. 

Next we define the total SNR for one triangle unit by 



0=1,11-'/"*. bh(I) 

and it has the following relation with respect to the geometrical parameters 



(28) 



SNR Z cx 



1 + cos i 



cos 2 i 



1 + cos 2 6 S 



cos 2 (2'0s) + cos 2 S sin 2 (2V>s) 



1 + cos 2 S 



sin 2 (2-0 s ) + cos 2 S cos 2 (2V>s) 



(29) 
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where i is the inclination angle defined as cos i = L-N . It is important to note that the right-hand side of this relation 
does not depend on the azimuthal angle cf>s- This is related to the orthonormality of the two effective detectors I and 
II. 

Up to 2PN order, we obtain the relation between the GW frequency / and the time t as 



t(f) = t c -^M z (KMJ)- 8 / 3 



, 4 /743 11 \ 32 , /2 

1 + 1 X) X 7TX 

3 V 336 4 7 5 

/ v -, 30 

„ / 3058673 5429 617 2 \ 2 1 v ' 



V 1016064 1008 144 

Here t c is the coalescence time, 77 = mxmz/imi + tti 2 ) 2 is the reduced mass ratio, and x is defined as x = (7r(l + 
z)M f) 2 / 3 with the total mass M = mi + 777,2- Since we neglect NS spins in this paper, we have omitted the spin- 
orbit coupling and the spin-spin coupling. When we only take the leading part and set t c — 0, the time before the 
coalescence is given by 

For z = 0, 3 and 5, the time —t(f — 0.2Hz) becomes 3.32 x 10 7 sec, 3.29 x 10 6 sec and 1.67 x 10 6 sec, respectively. In 
Fig. [5j we plot the accumulated averaged SNRs (normalized appropriately) against the time before coalescence for 
NS/NS binaries at z = 0,1,3 and 5. Since there would be WD/WD binary confusion noises below / = 0.2Hz (see 
Fig. [3]), we performed calculations from / = 0.2Hz to / = 100Hz. From this figure, it can be seen that the effective 
observation time depends strongly on the source redshift. There is almost 1 yr observation time for a binary at z = 0, 
but it is less than a month for a binary at z = 5. 

In reality, each unit of DECIGO or BBO is moving around the sun with the period T = lyr and characterized by 
the position angles 

0(t) = tt/2, 4>{t) = 2-Kt/T + c (32) 

with a constant Co (in this paper we take Co = 0). As we described in the previous paragraph, the observation times 
become longer for lower redshift binaries and it is important to take the effect of detector motions into account to 
perform more practical analyses. The time dependent amplitude v4 po i iQ essentially corresponds to the right hand side 
of Eq. (|2T)|) which is a function of S , ip s and i in the moving detector frame. We can change these variables to the 
angles (#Sj ^Si^Lj^l) i n the fixed frame by using the formulae shown in Appendix. 

Next we summarize the angular averaged SNR for binaries at a redshift z observed by a network of totally iV un ;t 
units of detectors, assuming that each contains two independent interferometers. Using the simple identity (a 2 o1 a 
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FIG. 5: This figure shows the accumulated averaged SNR of each (1.4+1.4)M© NS/NS binary against time to coalesce. The 
(red) circular plots, the (blue) triangular plots, the (green) square plots and the (black) crosses correspond to the one with the 
source redshift z = 5,3,1 and 0, respectively. These values have been normalized so that the accumulated SNRs become 1 at 
the time corresponding to /=100Hz. 



(4/5) 2 for the angular average, we have 

/•/max f-7/3 

(SNR 2 Z ) = 8N unit A 2 j -g^jdf- 

For a given position, the SNR of a binary becomes minimum for an edge-on geometry, and we have ^^4p i a 
1/5 for this subclass. Then the averaged SNR for these binaries is given as 



(33) 



edge— on 



pi ee (SNRf) , = — (SNRl) 

' z \ z l edge— on lfi 



(34) 



and is independent on the network geometry. Hereafter we mainly characterize the total sensitivity of a detector 
network by the parameter denned at the most distant redshift z = 5 in our fiducial model. This would be a 
convenient measure for discussing a network sensitivity required to identify foreground binaries. With iVunit = 4, 
/min = 0.2Hz and / max = 100Hz, we have the network sensitivity p$ = 10.4 for DECIGO and ps = 33.6 for BBO. 
The result p^ for BBO here is ~ 20% smaller than that in Cutler & Harms 
frequency cut-off / m i n - 

Since Eq. (|3"4")l is a mean value, a given network must satisfy a constraint 



14j . due to the difference in the lower 



SNR Z 

Ctz — 1 



(35) 



with respect to its minimum value SNR z jnin for binaries at a redshift z. The equality holds only when the SNR is 
independent on the direction and the polarization angles, e.g. for infinite number of units with random orientations 
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FIG. 6: (Left) The dead angle for gravitational wave interferometers with one detector plane. N is the unit vector from the 
detector to the binary, the unit vector z is normal to the detector plane and the unit vector L represents the orientation of 
the orbital angular momentum of the binary. L is perpendicular to N and has an angle 45° from z. We also introduce a new 
Cartesian coordinate (a/, y' , z') at the detector with x' = N and z' — z. (Right) The linear polarization pattern on the y' — z 
plane for the GW signal from a dead-angled binary (note that the other orthogonal polarization mode vanishes for the edge-on 
binary). An arbitrary detector arm in the x' — y plane is completely blind to this signal. 



(see also 3l|). Note that the ratio G z does not depend on the overall sensitivity of the detector (but on the shape of 
its noise spectrum). 



ANALYTIC STUDY FOR THE MINIMUM SNR 



Our primary interest in this paper is whether DECIGO or BBO can successfully identify the cosmological NS/NS 
binaries in performing their subtraction with sufficient accuracy to uncover an underlying stochastic background. The 
critical aspect here is the profile of the lower end of the probability distribution function of their SNRs. While we 
extensively perform numerical analyses in the next section, the simple analytical studies in this section would help 
us in understanding basic geometrical relations and also to interpret our numerical results. In this section, we fix the 
redshift z of binaries and do not include the annual motions of detectors for simplicity. The latter prescription would 
be a reasonable approximation for binaries at a large redshift z. 

Under these conditions, we discuss the minimum value SNR zul i n mostly dealing with edge-on binaries. As we see 
below, we have SNR Zl min = with one detector plane, due to the apparent dead angle described in Fig. |6l With 
multiple detector planes, the dead angle generally disappears, resulting in SNR Zjm i n > 0. Thus we evaluate the 
performance of a network with multiple planes, using the normalized measure G defined in Eq. (|35p (in this section, 
we simply denote G z as G since it does not depend on z for the static detector networks). The basic question here 
is how far we can increase the value G with a small number of units. As concrete examples, we pick up the specific 
configurations shown in Fig. [7] 
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FIG. 7: Configurations made by four units of detectors. In reality, each triangle unit is inclined by 60° from the orbital 
plane. (Top Left) We put all units on one site. There exists only one detector plane for this configuration. We label this as 
configuration "CI". (Top Right) We place three units on one site and put the remaining one with opening angle 02- There are 
2 detector planes for this configuration "C2(i)" . (Bottom Left) We place two out of four units on one site and the remaining two 
units on another site with a separation angle of 02- We label this as configuration "C2(ii)". (Bottom Right) This configuration 
"C3" is the same as the one shown in Fig[T] There exist three detector planes for this configuration. 

A. one detector plane 

As a first step, let us consider the situation only with a single detector plane (CI in Fig. [7]). From Eq. ([2T)|). we 
can straightforwardly confirm that the total SNR becomes identical to zero, when the following 3 conditions are 
simultaneously satisfied; cosi = 0, #s — 7r /2 and cos2^g = 0. This configuration is illustrated in the left panel of 
Fig. [6j Here, we introduce a new Cartesian coordinate (x',y',z') with x' = N and z' = z. For an edge-on binary, 
there is only one linear polarization mode (see Eqs. (|25|) and (|27j) ). Therefore we focus on the relevant mode shown 
in the right panel of Fig. [6] On the y' — z' plane its principal axis is tilted by 45° from the normal vector z' . We can 
evaluate the response of a detector arm with an angle 6 aim from the direction x' as [21 1 










f COS #arm\ 


(■ 





■ 


sin6> arm 




1 




^ ; 



0. (36) 

Since this equation holds for an arbitrary # arm , any GW detector on this plane is completely blind to this incoming 
signal. 

By expanding Eq. ([29]) around the dead angles (9s,ips,i) — ( 7r /2, 7r/4, 7r/2), we get, up to the leading order, 

9 dOl d%bl di 2 , , 

p 2 cx — £ + + — 37 
4 4 4 

with (d&s, dips,di) = (9s — tt/2, ■0s — 7r/4, i — ir/2). We will revisit this simple model later in Sec. IVI1 
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FIG. 8: (Left) ai denotes the separation angle between two detectors and fi% represents the angle formed by the planes on 
which these two detectors lie. (Right) Positions of two detectors and unit normal vectors z\ and Z2 are shown. Small panel on 
top left shows the y — z plane, indicating that each detector tilts 60° inwards from z axis. From the product z\ ■ £2, we have 
the correspondence cos {5% — (3 cos a.2 + l)/4. 



B. two detector planes 



Next, we consider networks with two detector planes, namely configurations C2(i) and C2(ii) shown in Fig. [7] In 
the present geometrical analysis, the network characterized by the position angle ai in the top left panel of Fig. [8] is 
equivalent to the network denned by the opening angle fa given in the bottom left panel of the same figure. These 
two angles are related with 

3 COS tt2 + 1 



COS ,62 = 



(38) 



which is derived by considering the inner product of the unit normal vectors Z\ and £2 for each detector plane shown 
in the right panel of Fig. [3] 

With the dead angles for the previous subsection, it can be easily understood that we have G = only for the 
specific angles fa — (mod tt/2). For other angles, we find that the minimum value SNR Z:m i n (> 0) is given by an 
edge-on binary at the direction of the intersection of two planes. From Eq. (|2U)) we can show that, as a function of 
fa, the maximum value of G is realized at the optimal angle fa = tt/4 with 



G(fa = tt/4) 



2^2' 



(39) 
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FIG. 9: The normalized ratio G = SNR z ,uiin/pz for C2(i) and C2(ii) as functions of the opening angle /?2- The optimal angle 
is /?2 = tt/4 (corresponding to tt2 = 52.4°). 

Here the parameter r(< 1/2) represents the relative weight of the number of detector units between two sites (e.g. 
r = 1/4 for C2(i) and 1/2 for C2(ii)). In Fig. [9] we show the ratio G for the configurations C2(i) and C2(ii). The 
optimal angle f3 2 = 7r/4 corresponds to a 2 = 0.915 (52.4°), and we have G(-k/4) = V5/4\/2 = 0.395 for C2(i) and 
a/5/4 = 0.559 for C2(ii). For the latter C2(ii), we can further derive a simple expression as 

Gifa) = 0= min[| cos fo\ , | sin |] . (40) 

As one can see, when a 2 equals to tt (maximum separation), the opening angle (3 2 becomes 27r/3 and this configuration 
is not the optimal choice. 

Although Fig. [9] is symmetric with respect to the value f3 2 — 7r/4, this does not mean that the situations are the 
same for f3 2 < ir/4 and (3 2 > ir/4. For example, let us focus on two configurations with opening angles of (3 2 = 
and 7r/2. j3 2 — corresponds to a configuration with only one detector plane while there are two detector planes for 
f3 2 = it/2. Since each has some binaries that are completely insensitive to each configuration, they both give G = 0. 
However, the fraction of these dead-angle binaries to whole binaries is smaller for /3 2 — tt/2 than the one for f3 2 = 0. 
Therefore we expect that the symmetry of Fig. [9] is broken when we take the motions of detectors into account. 

C. three detector planes 

As an extension of the default configuration sketched in Fig. [TJ we relax the constraint 03 = 120° for the angle 
between the three positions of units, and evaluate the normalized ratio G as a function of a^. Our result is presented 
in Fig. [TO] The ratio G takes its maximum value 0.771 at the optimal angle 0:3 = 2.205 (126.36°) that is slightly 
larger than the default angle 0J3 = 120°. 
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FIG. 10: The normalized ratio G = SNR Zim x n /p z for C3 as a function of the detector separation angle 03 . We have G = 0.771 
at the optimal angle 03 = 2.205 (126.36°). 

For comparison, we examine a similar network but with evenly weighted sensitivity (namely 1:1:1 instead of 2:1:1 
as in Fig. QJ. In this case, the maximum of the ratio G becomes 0.843 at the symmetric configuration 013 = 120°. It 
seems reasonable that, as we increase the weight of the overlapped units shown in Fig. [TJ the optimal angle shifts 
to a larger value from the symmetric value 013 = 120°. We also evaluate the parameter G when more than three 
detectors are available on the ecliptic plane and evenly weighted. In the large number limit, the asymptotic value 
can be given analytically as G = V205/16 = 0.895 



231 ] . Note that even if we include the detector motions, this value 
is unchanged and thus can be regarded as an useful upper bound for the BBO/DECIGO-type orbital configurations. 
We expect that as we increase the symmetry of a detector configuration, the effect of detector motion is suppressed. 

One might be interested in a network formed by evenly weighted three units with each normal to the three orthogonal 
directions (1,0,0), (0,1,0) and (0,0,1). Even though we do not have a direct orbital model for space detectors, it 



could be an efficient network in terms of the ratio G. Indeed, we have G = 5/v30 = 0.913 for this hypothetical 



network (see also [3l|). 

In summary, using four equivalent units shown in Fig. [T]with a free parameter 03, we can realize G = 0.771 that is 
at least 39% larger than a network with two detector planes. The optimal angle CV3 is 126.36° and somewhat larger 
than the default value 120°. Here, we should remind that our results in this section are derived by neglecting motions 
of detectors, and valid for the binary sources at the high redshift limit. In the next section, we numerically evaluate 
the parameter G for realistic cosmological binaries. These binaries have finite duration times in the BBO/DECIGO 
band (see Fig. [5} and we should include the detector motions when studying their SNRs. 
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VI. NUMERICAL CALCULATIONS AND RESULTS 
A. Set Up 

In order to estimate how well we can subtract the foreground GWs by NS/NS binaries, we performed Monte 
Carlo simulations in the following manner. First we uniformly divide redshift range z = [0, 5] into 200 bins. For each 
redshift bin, we randomly generate 10 sets of binary directions (9s,4>s) and orientations (#l,</>l)> with cos 9s and 
cos^l uniformly distributed in the range [—1,1] and </>s and 0l in the range [0, 2ir}. Then we calculate SNR of each 
binary, fixing the redshifted NS masses at 1.4(1 + z)M®. We denote the SNR of i-th binary out of the 10 4 binaries at 
redshift z as p(z, i). We start integrations of GW signals from / m j n = 0.2Hz, considering the potential confusion noise 
by WD/WD binaries. We set the upper cutoff frequency / max = 100Hz, but our results are insensitive to this choice. 

For the detector configurations, we consider 4 cases shown in Fig.[7]with DECIGO. To concentrate on the geometrical 
effects, we use the fixed number of the triangle units (totally four, corresponding to eight L-shaped interferometers). 
These units are assumed to have the identical sensitivity to DECIGO, and only their positions and orientations are 
different. Note that our analysis is based on the long- wave approximation for responses of detectors, and the two 
units forming a star-of-David constellation can be regarded as a completely aligned pair. The results below remain 
almost the same when we use BBO instead (as we confirm in Sec. IVICp . 

B. Minimum SNR 

In Fig. [TTJ we show the probability distributions of SNRs for the edge-on NS/NS binaries at z = and 5. We 
count the number of binaries within each SNR bin and divide it by the total number of binaries 10 4 . In order to 
compare the overall distribution pattern of SNRs at the two redshifts, wc normalize the SNRs by p z . We consider 
the configurations CI, C2(i) and C3, with the latter two at the optimal separation angles {a^ — 52.4°, — 126.36°). 
For reference, we also consider a somewhat extreme configuration denoted as CI'. This has totally eight L-shaped 
interferometers placed on one site. 

First, let us look at the results of z — 5. We see that as we increase the number of detectors, the distributions 
become sharper due to the averaging effect of the beam pattern functions. For the distributions of CI and CI', there 
are tails on lower SNR side, down to SNR^/p^ m 0. This is because there exist dead angles of binaries for the static 
detectors having only one detector plane. However, since the detectors are moving during the observation period of 
about 3 weeks, the minimum SNR for each configuration is not exactly 0. On the other hand, there exists a distinct 
cutoff on lower SNR side (and also on higher SNR side) for the configurations C2(i) and C3. This is because these 
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FIG. 11: These figures show the probability distributions of binaries in each SNR bins, comparing the results for both binaries 
at z = and z = 5. We normalize each SNR with the root-mean-square value p z . Each panel on top left, top right, bottom 
left and bottom right corresponds to the one using detectors of the configurations CI', CI, C2(i) and C3, respectively. Here, 
CI' represents the configuration with 8 right-angled interferometers placed on one site, as shown in the small panel on top left. 
The (red) triangular plots in each panel show the results for binaries at z = 5 while the (blue) crosses indicate the ones at 
z = 0. The solid line in each panel of C2(i) and C3 represents the value G for each optimal configuration with a static detector 
network. The dashed line in C3 panel shows the one G = 0.895 with infinite number of detectors (see Sec. IV Cj) . 



networks with more than 1 detector plane do not have dead angle directions even for the static case. This cutoff 
represents the minimum SNR (normalized with p 5 ) for all NS/NS binaries detected by DECIGO/BBO and therefore 
corresponds to the value G in Sec. [V]but now for the case of moving detector networks. If the detector sensitivity is 
high enough to identify the signal of this lowest SNR binary, we can handle individual contributions of all binaries 
in the data streams of detectors (provided that, ultimately, subtraction of individual binaries are not limited by the 
NS/NS foreground itself and the residual noises after the binary subtraction would be sufficiently small). In each 
panel of C2(i) and C3 in Fig. [Til we show the analytical value of G for the static detector network as a solid line, 
which is obtained in Sec.[V] We can see that this analytical value almost matches with the lower cutoff value obtained 
by the Monte Carlo simulation. Again, due to the motions of detectors, the numerical counterparts are slightly larger 
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FIG. 12: The normalized ratio G5 = S TV i?5, m in / 'ps for C2(i) and C2(ii) as functions of the opening angle 02- The solid line 
represents the analytical results shown in Fig. [9] assuming that detectors are static. The dots are the ones obtained by our 
numerical calculations which include the effects of detector motions. 

than the ones estimated analytically for the static configurations. 

Next, we compare the results of z — 5 and z — 0. Figure [Til shows that the SNR probability distributions of z = 
are sharper than the ones of z = 5. This is because the observation time for z = is longer than for z — 5, and 
the number of independent detector planes arc effectively increased. In the panel of configuration C3, we also show 
the analytic value G = v / 205/16 = 0.895 given for infinite number of detectors as a dashed line. This value almost 
matches with the lower cutoff value of z — in the same panel. Since observation time of binaries at z — is slightly 
larger than 1 yr, detectors orbit around the Sun completely. Therefore, as we increase the number of detector planes, 
the lower cutoff in the probability distribution approaches to G — 0.895. 

Now, we compare the analytical results of the values G for the static detectors shown in Figs. [5] and [TU] with our 
numerical results which include the effects of detector motions. In Fig. [T^l we present G5 for the configurations C2(i) 
and C2(ii) that have only two detector planes. Solid line on each panel indicates the analytical result, while the dots 
are the ones obtained by our numerical simulations. We can see that the analytical and numerical results are quite 
close for the cases of large G5, while they do not match for the cases of small G5, especially where G5 « 0. This shows 
that the effects of detector motions are more significant for the smaller G5 cases. When detectors are static, there 
exist dead angles for the configurations /?2 = and /?2 = 7r/2. These dead angles disappear when detectors are moving, 
which makes G5 =/= and causes a remarkable difference between the static and the moving cases. Also, numerical 
results are not symmetric with respect to P2 — 7r/4, as discussed in Sec. V.B. For the case of $2 = 0, the number 
of detector plane becomes one and the results for C2(i) and C2(ii) should match. The discrepancy comes from the 
statistical fluctuations which are caused by the finiteness of the number of binaries in our Monte Carlo simulations. 
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FIG. 13: The normalized ratio G5 = SNRs tU1 i n /p5 for C3 as a function of the detector separation angle Q3. The meanings for 
the line and dots are the same as in Fig. 1121 

In Fig. 1131 we show the results for the case C3. In this case, the distinction between analytical and numerical results 
is smaller than the cases C2. The maximum value of G5 is obtained when 013 is around 126°, though it is almost the 
same with the one for the default configuration of 03=120°. From Figs. [T2l and [T3l we understand that the matching 
between the analytical and numerical results becomes better as we increase the symmetry of detector configurations. 
This tendency is consistent with the discussion made in Sec. [V] As we increase the symmetry, the effect of detector 
motions is suppressed and in the limit of infinite number of detector configuration, the values of G5 are exactly the 
same for the static and moving cases. 

Next we evaluate the network sensitivity P5 required for identifying all the binaries with our fiducial merger model 
in which NS/NS signals exist up to z — 5. With respect to the minimum value SNRs^ m - m of our binaries, we have 

SNR 5 

.min 



P5 



(41) 



due to the definition of the parameter G5. For complete identification of the binaries, we should have 



SNR 5:min > p t hr 



(42) 



with the detection threshold pthr (which 



required to cover the parameter space 



may depend on the future computing power and the number of templates 



14J). Then we get 



Pth 

P~5 



L <G 5 



(43) 
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FIG. 14: (Left) This figure shows the accumulated fractional residuals R against pthr/ps for various configurations. The (red) 
circular plots, the (blue) crosses and the (black) triangular plots show the ones using detectors of configuration CI, C2(i) and 
C3, respectively. For each case of C2(i) and C3, we take the optimal configuration. We assume that the subtraction is succeeded 
if the residual is below the horizontal dashed line of R = 10 . (Right) This figure also shows the residuals R for the detectors 
having only 1 detector plane. The (red) circular plots are the results using the moving DECIGO and this is same as the one 
shown in the left panel. The (blue) crosses indicate the ones using the moving BBO while the (green) square plots are the ones 
using the static DECIGO, with the fitted line of R oc (pthr/ps) shown as the black solid line. 



or equivalcntly 



P5 > G 5 Vthr = 25.9 



G 5 



( fkhr 



(44) 



0.77 ) V 20 

where we used the parameter G5 = 0.77 for the case C3 and the typical threshold pthr = 20. With the network 
sensitivity of p<- = 10.4 for DECIGO, many binaries have SNRs below pthr = 20 and would not be identified. To 
fix this, the sensitivity of DECIGO should be improved by a factor of 2.5 (^fj. In contrast, BBO has ps = 33.6 
and, with G 5 = 0.77, it can detect all the NS/NS signals with SNRs higher than 33.6x0.77=25.9. However, for a 
higher threshold pthr = 30, the network sensitivity is not sufficient even for BBO and then we are not able to perform 
complete subtraction. In the next subsection, we study what fraction of them remains to be unidentified. Notice that 
the results in this subsection are not affected by the values of no- 



Residuals 



In order to estimate the amplitude of unidentified foreground composed by binaries with SNRs less than the 
threshold pthr, we define the ratio 7(2) for the numerical samples described in Sec. VLB as 

£iP(M) 20 [Ahr ~ P{z,i)} 



7(2) 



EiP(M) a 



(45) 
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with the step function ©(•). Then, from Eq. ((8|), we can calculate the fraction of the residual NS/NS binaries in the 
total foreground as 

Jo d»7W (1 y gw 

r°° r/ y s ( z ) 1 
JO a<i (l+ 2 ) 4 / 3 H(^) 

The basic question here is how the residual fraction R depends on the configuration of detectors. 

In the left panel of Fig. [TH we show the residual R against pthr/ps for the cases CI, C2(i) and C3. For the case of 
C3, the residual falls down to very steeply as we lower pthr/f>5- This figure shows that in order to achieve R < 10~ 5 
for our fiducial merger rate of no = 10~ 7 Mpc~ 3 yr _1 , the network sensitivity ps should satisfy pthr/ps < 0.91. 
However, since the value of R depends very sensitively on pthr/P5j we suggest that it seems much safer to aim for 
Pthr/f>5 < 0.77 so that all of the binary signals can be identified. For the case of C2(i), the residual R depends more 
gently on pthr/ps and R < 10 requires pthr/p5 < 0.63. In this case, it might be reasonable to adjust the network 
sensitivity p$ to perform successful subtraction down to R < 10 . However with the sensitivity pthr/ps < 0.63, we 
suggest to use the C3 configuration and subtract all the binary signals. For the case of CI, we need a better sensitivity 
Pthi/p5 < 0.48 to accomplish R < 10~ 5 , which seems a rather demanding task. 

Next, we compare the results for DECIGO and BBO to check that they give almost the same results. Out of the 
configurations CI, C2 and C3, it is expected that the configuration CI gives the largest discrepancy between these 
two results. On the right panel of Fig. [14J we show the dependence of residuals on pthi/f>5 using detectors having 
only one detector plane (CI). The (red) circular plots are the ones using DECIGO and this is the same as the ones 
on the left panel, while the (blue) crosses represents the ones using BBO. Since the results for the moving DECIGO 
and BBO are almost the same, we can confirm that our main results obtained in this paper can be applied to BBO 
as well. 

On the same panel, we show the residuals obtained by using the static DECIGO and also provide a fitted line of 
R oc (pthr/ps) 5 - This index of 5 can be understood analytically, as follows. When considering detectors with one 
detector plane, we have the dead angles (9s, ips, i) = (f/2) 7r /4, f/2) and we expect that the unidentified binaries have 
angular parameters similar to these specific combinations. Thus, around the dead angles, we introduce the coordinate 
system (A0g, Alps, Ai), and consider the region where the SNR is less than a given threshold pthr- From Eq. (|37[) . 
the region turns out to be a sphere, and its radius is proportional to pthr- Therefore, the number of binaries -/V <t hr 
within the region has the following relation 

7V<thr oc A0 s A^ s Ai cx p 3 thr . (47) 
Since the averaged SNR in the region is proportional to pthr, together with Eq. (|4"7|) . the unidentified fraction defined 
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in Eq. (l45l scales as 

7 oc p t 5 hr . (48) 

Since the fractional residuals R is roughly proportional to 7, we can reproduce the simple dependence R oc (pthr/ps) 5 - 
When we include the detector motions, the residual starts to deviate from this power-law as we lower the value of 

Pthr/p~5- 

VII. CONCLUSIONS 

In this paper, we estimated the detector noise level required to individually identify foreground GWs made by 
cosmological NS/NS binaries. We characterized the total sensitivity of the detector network by an integral ps, 
assuming that the binaries exist at z < 5. We mainly focused on the relationship between the geometrical properties 
of the detectors and the minimum SNR of the binaries, and introduced a convenient parameter G z defined by 

G z = SNR *<™ < 1. (49) 

Pz 

We started our studies with static detector networks. In contrast to a network with a single detector plane, in general 
the minimum SNR of binaries becomes finite for the one with more than 1 detector planes. For the C3 configuration 
(see Fig. [1} with 3 detector planes, we analytically estimated the optimal angle 0:3 = 126.4° with G z = 0.771. It is 
the angles between detector planes that play important roles in realizing a high G and not the separations between 
detector units (although they are important for improving the angular resolution). 

We expected that the ratio G z becomes larger when we use the moving detectors, due to an averaging effect. In 
order to estimate this correction, we performed Monte Carlo simulations, randomly distributing cosmological NS/NS 
binaries at z < 5. We showed that our analytical predictions and numerical results match well, especially for large G5 
cases. In order to identify all of NS/NS binaries, we need to achieve enough detector sensitivity so that the minimum 
SNR is larger than the detection threshold pthr- This leads to the following inequality 

«> G r'"«„ = 25.9( j ^)" 1 (^). (50) 

It suggests that for a given threshold at p t hr = 20, all the binary signals can be identified if the optimal C3 network 
has a sensitivity p~5 > 25.9. With its reference parameters, DECIGO has the network sensitivity ps = 10.4. In order to 
perform complete identification of NS/NS binaries at z < 5, its sensitivity should be improved by a factor of 2.5 (^-) • 
BBO has ps = 33.6 and can identify all of them as long as pthr = 20. However, for a higher threshold pthr = 30, even 
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BBO fails to subtract all the binary signals. At the end of our calculation, we estimated the residual fraction R of 
the foreground that cannot be identified. We found that the low SNR end of the residual is very steep. Therefore we 
recommended to realize the sensitivity sufficient for complete identification rather than those with R < 10~ J . 

When calculating SNR in this paper, we only used the instrumental noise spectral density. In order to perform more 
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realistic analysis, we need to include NS/NS foreground noises and perform self-consistent analysis as done in Ref. 
Furthermore, although the individual GW signals are assumed to be removed once they have been identified, this is 
not a straightforward step. With respect to this issue, we should perform more elaborate data analysis simulations 



(see e.g. Ref. 



15j and also Refs. 
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14| discussed that typical NS/NS 



In this paper, we have assumed circular binaries with no spins. Cutler and Harms 
binaries have squared eccentricities e 2 ~ 10~ 8 at / = 0.3Hz. For these binaries, the correction in the GW phase due 
to eccentricities become ~ 0.2 for the dominant harmonic, and we should include this phase effect in the template for 
matched filtering. They also showed that the spin-orbit correction to the phase becomes important when one of the 
binary component has spin period ~10 msec. However, even in this case, the contributions from spin-spin couplings 
and precessions can be neglected. 

It is expected that we can strongly constrain cosmological parameters by observing GWs from chirping binaries. 
The crucial point here is whether we can determine the redshifts z of the GW sources by associated electro-magnetic 
wave observations (see Cutler and Holz [24| and references therein). It might be useful to study dependence of the 
angular resolution of binaries on configurations of detectors. Another aspect closely related to the geometry of the 
detector network is the decomposition of various polarization modes of GWs, including parity asymmetry (the Stokes 



V parameter) (34j and the non-Einstein modes 35 [. 
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Appendix A: Useful formulae 

When we perform Monte Carlo simulations, we randomly distribute the direction of the source (9s,<j>s) and the 
direction of the orbital angular momentum (#l> </>l), both measured in the solar barycentric frame. Therefore we need 
to express the waveforms in terms of 6 S , S , # L and 0l (see [21J for details). The angles 6>g(t) and 4>s(t) m the moving 
detector frame are expressed as 

1 - V3 

cos6» s (i) = -cos6» s —sm9 s cos[4>{t)-<f>s], (Al) 

Mt) = JL + tan -i f ^cosg s + sin^cos[#) -0 s ]\ (A2) 
VW 12 ^ 2sin(9ssin[0(t)-0 s ] J y ' 

The polarization angle "0s in the moving detector frame is given by Eq. (|24p with the following three quantities L ■ z, 
L ■ N and N ■ (L x z). Here we have 

1 - V3 

L ■ z = - cos 6» L r~ sin 6> L cos [</>(*) - L ] , (A3) 

L-N = cos^l c°s^s + sin0 L sin^s c °s(^L — ^s)) (A4) 
N-(Lxz) = i sin L sin S sin(0 L - S ) 

— — cos</)(t)(cos0L sin0 s sin^g — cos^s sin^L sin<^ L ) 

— — sin0(t)(cos^s sin^L cos^l — cos#l sin 9s cos^s). (A5) 
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